Watersheds are Schramm-Loewner Evolution curves 
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We show that in the continuum limit watersheds dividing drainage basins are Schramm-Loewner 
Evolution (SLE) curves, being described by one single parameter k. Several numerical evaluations 
are applied to ascertain this. All calculations are consistent with SLE K , with k = 1.734 ±0.005, 
being the only known physical example of an SLE with k < 2. This lies outside the well-known du- 
ality conjecture, bringing up new questions regarding the existence and reversibility of dual models. 
Furthermore it constitutes a strong indication for conformal invariance in random landscapes and 
suggests that watersheds likely correspond to a logarithmic Conformal Field Theory (CFT) with 
central charge c w — 7/2. 
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The possibility of statistically describing the proper- 
ties of random curves with a single parameter fascinates 
physicists and mathematicians alike. This capability is 
provided by the theory of Schramm-Loewner Evolution 
(SLE), where random curves can be generated from a 
Brownian motion with diffusivity ft Once re is iden- 
tified, several geometrical properties of the curve are 
known (e.g. fractal dimension, winding angle, and left- 
passage probability) [21 [3] . Among the examples of such 
curves, we find self-avoiding walks [4j and the contours 
of critical clusters in percolation [5] , Q-state Potts model 
[5], and spin glasses [7], as well as in turbulence [SJ. Es- 
tablishing SLE for such systems has provided valuable 
information on the underlying symmetries and paved the 
way to some exact results [51 HUTU]. In fact, SLE is not 
a general property of non-self-crossing walks since many 
curves have been shown not to be SLE as, for example, 
the interface of solid-on-solid models [TT], the domain 
walls of bimodal spin glasses [13], and the contours of 
negative- weight percolation (13j . 

Recently, the watershed (WS) of random landscapes 
[14TfT6] . with a fractal dimension df ~ 1.22, was shown to 
be related to a family of curves appearing in different con- 
texts such as, e.g., polymers in strongly disordered media 
[T7j , bridge percolation jT3] , and optimal path cracks [TH] • 
In the present Letter, we show that this universal curve 
has the properties of SLE, with re = 1.734 ±0.005. re < 2 
is a special limit since, up to now, all known examples 
of SLE found in Nature and statistical physics models 
have 2 < re < 8, corresponding to fractal dimensions df 
between 1.25 and 2. 

Scale invariance and, consequently, the appearance of 
fractal dimensions have always motivated to apply con- 
cepts from conformal invariance to shed light on criti- 
cal systems. Archetypes of self-similarity are the con- 
tours of critical clusters in lattice models. Already back 
in 1923, Loewner proposed an expression for the evo- 
lution of an analytic function which conformally maps 



the region bounded by these curves into a standard do- 
main |19j . Such an evolution, follows the theory, should 
only depend on a continuous function of a real parame- 
ter, known as driving function. Recently, Schramm ar- 
gued that to guarantee conformal invariance, and do- 
main Markov property, the continuous function needs 
to be a one-dimensional Brownian motion [T] thrusting 
into motion numerous studies in what is today known 
as the Schramm-Loewner, or Stochastic-Loewner, Evolu- 
tion (SLE). Such one-dimensional Brownian motion has 
zero mean value and is solely characterized by its diffu- 
sivity re, which relates with the fractal dimension df as 

EDI EH, 

d f = min{l + re/8, 2} . (1) 

Although it is believed that SLE should hold for the en- 
tire class of equilibrium O(n) systems, it has only been 
rigorously proven for a few cases [SJHUj. Nevertheless, 
numerically correspondence has been shown for a large 
number of models as mentioned above. It has been ar- 
gued that SLE can be applied to models exhibiting non- 
self-crossing paths on a lattice, showing self-similarity, 
not only in equilibrium but also out of equilibrium as the 
example discussed here [SJ [33J . 

In random discretized landscapes each site is charac- 
terized by a real number such as, e.g., the height in an 
elevation map, the intensity in a pixelated image, or the 
energy in an energy landscape [2]. If sites are occupied 
from the lowest to the highest, clusters of adjacent oc- 
cupied sites can be defined and, at a certain fraction of 
occupied sites, a spanning cluster emerges connecting op- 
posite borders (e.g. from left to right). In the example 
of the elevation map, this procedure corresponds to fill- 
ing the landscape with water until a giant lake emerges at 
the threshold, which drains to the borders [24] , When we 
now suppress spanning by imposing the constraint that 
sites, merging two clusters touching the opposite borders, 
are never occupied, a line emerges delineating the bound- 
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FIG. 1. (color online) Denning on as the turning angle 
between two adjacent edges, i and i + 1, their winding angles 
8i and 9i+i are related by 9i+\ = 9% + oti, with 9i = 0, as 
illustrated in the main figure. We performed simulations for 
18 different lattice sizes, with L y = 2 4+n and 3 x 2 3+n for 
n — 1,2, ... ,9. Results are averages over 10 6 samples for the 
smallest system sizes and 3 x 10 3 for the largest one. Inset: 
Dependence of the variance of the winding angle on the lateral 
size of the lattice L y . Statistical error bars are smaller than 
the symbols. The slope in the linear-log plot corresponds to 
K /4 = 0.432 ± 0.002. 

aries between two clusters: one connected to the left and 
the other to the right border [14]. This line is the wa- 
tershed line (WS) separating two hydrological basins [25] 
and has a fractal dimension df = 1.2168±0.0005 [2B]. We 
show in this Letter that in the scaling limit its statistics 
converges to an SLE K , consistent with k = 1.734± 0.005. 

To study the scaling limit of WS and compare it to 
SLE K , we numerically generated ensembles of curves and 
carried out three different statistical evaluations, namely, 
the variance of the winding angle (quantifying the angu- 
lar distribution of the curves) [2ZJ [25], the left-passage 
probability [TJ[25], and the characterization of the driving 
function (direct SLE) [8j. We show here that the values 
of k independently obtained for each analysis are numer- 
ically consistent and in line with the fractal dimension of 
the WS. For all cases, simulations have been performed 
on both, square lattices of square shape (L x = L y ) and in 
strip geometries (L x > L y ), all with free boundary con- 
ditions in horizontal and periodic boundary conditions in 
vertical direction. L x is the size of the horizontal bound- 
ary, while L y is the length of the vertical one. Hereafter, 
we discuss each analysis separately. 

Winding angle. Using conformal invariance and 
Coulomb-gas techniques, Duplantier and Saleur [27] have 
found the dependence of the distribution of the winding 
angle on the system size and the Coulomb-gas parameter. 
Given the correspondence of the Coulomb-gas parameter 



FIG. 2. (color online) Original (left) and mapped (right) 
watershed. The Schwartz-Christoffel transformation (SCT) 
has been applied to map from the square lattice (left) to the 
upper half-plane (right). In this way dipolar curves are turned 
into chordal curves. 

to «, the relation for the winding angle can be extended 
to SLE [28]. To analyze the winding angle 6*.; at edge i, we 
set 9i = and define on as the turning angle between the 
edges i and i + 1 (see Fig. [l}. The winding angle of each 
edge is then computed iteratively as Qi+\ = 0i + on. For 
SLE K , the variance of the winding angle over all edges in 
the curve scales as (9 2 ) = b + (ft/4)lnL y , where & is a 
constant and L y the lateral size of the lattice [30] . Fig- 
ure [T] shows the variance as function of lateral size L y for 
the WS, with a slope 0.432 ± 0.002 in a linear-log plot. 
This slope corresponds to ft = 1.728 ± 0.008 which is in 
good agreement with the one predicted by Eq. ([I]) from 
the WS fractal dimension. 

From dipolar to chordal representation. In the original 
setup, WS are dipolar curves which start at one point 
on the lower boundary and end when they touch the up- 
per boundary, for the first time. For the left-passage 
probability and direct SLE evaluations, exact results are 
however known for chordal curves |29j , which start at the 
same point but go to infinity. Therefore, to proceed with 
these evaluations, we map the dipolar WS curves into 
chordal ones in the upper half-plane H (see Fig. |2|. For 
such mapping, as suggested in Refs. [3T1GI2], we usec ^ tne 
inverse Schwartz-Christoffel transformation |33j . 

Left-passage probability. For SLE curves in the upper 
half-plane H, starting at the origin, the probability that a 
point Re 1 ^ is at the right side of the curve (see Fig.[3|a)) 
solely depends on <f> and k and is given by Schramm's 
formula [29] . 

pm = cot ^^ Q; 2 ; - cot2 w ) • 
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FIG. 3. (color online) (a) Schematic representation of the 
left-passage definition (details in the text), (b) P((f>,R) — 
Pi. 725(0) for the chordal watershed at different distances from 
the origin R = {0.05,0.5,1.2}, where Pi. 725(0) is the left- 
passage probability for k — 1.725 given by Schramm's for- 
mula, Eq. j2|. (c) Mean square difference Q(k) between the 
numerical data and Schramm's formula (Eq. Q) for different 
values of k, exhibiting a minimum at K = 1.73 ±0.01. In both 
cases, results are averages over 10 5 curves on square lattices 
with L v — 512. 



where 2-F1 is the Gaussian hypergeometric function and T 
is the Gamma function. Figure [3jb) are the data points 
for the difference between the numerically measured 
probability P{<j>, R) and the one predicted by Schramm's 
formula, Eq. (J2j|, for the chordal curve. It is shown that 
P((f>,R) is independent on R. To estimate k we plot, in 
Fig. [3jc) , the mean square deviation Q(k) defined as, 



Q(k) 



(3) 



where the outer sum goes over values of 0.05 < R < 1.2, 
in steps of 0.05, and the inner one over values of < 
4> < 7r, in steps of 7t/15. M is the total number of con- 
sidered points Re % ^ . To reduce the statistical noise we 
used the relation P(4>, R) + P(ir — <j>,R) = 1. The mini- 
mum in the plot corresponds to the value of k that best 
fits the left-passage probability, giving k — 1.73 ± 0.01, in 
line with the prediction based on the fractal dimension 
of WS, given by Eq. Q. 

Direct SLE. Consider a chordal SLE curve j(t) which 
starts at a point on the real axis and grows to infinity in- 
side the region of the upper half-plane H, parametrized 
by an adimensional parameter t, typically called Loewner 
time. To compute its driving function one needs to 
find the sequence of maps gt(z) which at each time t 
map the upper half-plane EI into EI itself and satisfy the 
Loewner equation |19j . This map is unique and can be 
approximately obtained by considering the driving func- 



FIG. 4. (color online) Dependence of the second moment 
of the driving function (£ 2 (i)} on the Loewner time, for the 
chordal watershed. The slope corresponds to k = 1.69 ±0.05. 
Inset: Probability distribution of the driving function at two 
different Loewner times for chordal watersheds. The rescaled 
parameter X is defined as X — £(t)/\Act, where we have taken 
k — 1.69. The solid line is the normal distribution of vanishing 
mean value and unit dispersion. Results are averages over 
4 x 10 4 realizations on a square lattice with L y — 1024. 



tion to be constant within an interval 5t, obtaining the 
slit map, 



4^ 



(4) 



where z is a point in EI and 5t also depends on t. This 
map converges to the exact one for vanishing 5t [34] . Ini- 
tially, we set t = and £(0) = and we proceed iter- 
atively through all points of the chordal curve. At 
each iteration j, we map the point Zj to the real axis, 
by setting Stj — (Im Zj) 2 /4 and the driving function 
£(tj) — Re Zj (being Re and Im the real and imagi- 
nary parts, respectively). We also compute the Loewner 
time tj = tj—i + Stj. As referred above, in SLE, the 
driving function is related to a Brownian motion B(t), 
with vanishing mean value and unit dispersion, such that 
£(t) = y/RB(t) [35]. 

Figure [4] shows the second moment of the driving 
function for the chordal WS. The inset displays the 
probability distribution for the rescaled driving func- 
tion X = £(t)/y/Kt for two different times for the chordal 
WS. All results are consistent with a Brownian motion 
with vanishing mean value and unit dispersion, when 
K = 1.69 ± 0.05, in good agreement with the results dis- 
cussed above. The direct SLE analysis is characterized 
by larger error bars than the other two methods (wind- 
ing angle and left-passage probability) due to strong dis- 
cretization effects in the slit mapping [34] . 

Discussion. Our detailed numerical analysis shows 
that watersheds are likely to be SLE curves with 
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k = 1.734 ± 0.005. This is the first documented case of a 
physical model with k < 2, lying outside the well-known 
duality conjecture range 2 < k < 8, giving k' = 16/ k, 
where k' is the diffusivity of the dual model [35] . It has 
been proven that SLE K with k > 8 is not reversible [55] , 
therefore if a dual model exists which respects reversibil- 
ity, then it cannot be SLE K / with n' > 8. In the context 
of SLE, duality implies that two apparently different frac- 
tal dimensions might actually stem from the same curve. 
Geometrically, this corresponds to a relation between the 
fractal dimension of the accessible external perimeter and 
the one of the curve. 

Our work shows that watersheds are non-local SLE 
curves. Although a connection with SLE is strong indi- 
cation for conformal invariance, it cannot be interpreted 
as a proof. Nevertheless, if such invariance is established, 
it becomes possible to develop a field theory for this 
new universality class. CFT has helped to classify con- 
tinuous critical behavior in two-dimensional equilibrium 
phenomena 37, 38J . A well-established relation between 
diffusivity k and central charge c of minimal CFT mod- 
els which have a second level null vector in their Verma 
module is c = (3k — 8) (6 — k)/2k 34 . If the water- 
shed is conformally invariant it likely corresponds to a 
logarithmic CFT (LCFT) with central charge c « -7/2. 
A series of LCFT's corresponding to loop models have 
been suggested in Ref. [35], which thus seem to be re- 
lated to watersheds. It is also noteworthy that negative 
central charges have been reported in different contexts 
like, e.g., stochastic growth models, 2D turbulence, and 
quantum gravity 40-42]. In particular, the loop erased 
random walk is believed to have k = 2 which corresponds 
to c = —2. Besides, since the watershed of a landscape is 
based on the distribution of heights, the configurational 
space grows with TV!, where N is the number of sites, be- 
ing a promising candidate to develop a field theory with 
quenched disorder. 

The connection between SLE and statistical properties 
of the watershed opens up new possibilities. Since the 
latter are related to fractal curves emerging in several 
different contexts, our work paves the way to bridge be- 
tween connectivity in disordered media and optimization 
problems where the same k, and its corresponding central 
charge, are observed. Besides, a systematic study of the 
k dependence on correlations in the landscape might pro- 
vide the required information to find SLE curves on natu- 
ral landscapes. The possibility of a multifractal spectrum 
for watersheds is also an open question. 
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